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ABSTRACT 

We have performed magnetohydrodynamic (MHD) simulations of the collapse and fragmen- 
tation of molecular cloud cores using a new algorithm for MHD within the smoothed particle 
hydrodynamics (SPH) method, that enforces the zero magnetic divergence constraint. We find 
that the support provided by magnetic fields over thermal pressure alone has several important 
effects on fragmentation and the formation of binary and multiple systems, and on the prop- 
erties of massive circumstellar discs. The extra support suppresses the tendency of molecular 
cloud cores to fragment due to either initial density perturbations or disc fragmentation. Fur- 
thermore, unlike most previous studies, we find that magnetic pressure plays the dominant 
role in inhibiting fragmentation rather than magnetic tension or magnetic braking. In particu- 
lar, we find that if the magnetic field is aligned with the rotation axis of the molecular cloud 
core, the effects of the magnetic field on fragmentation and disc structure are almost entirely 
due to magnetic pressure, while if the rotation axis is initially perpendicular to the magnetic 
field, magnetic tension plays a greater role and can actually aid fragmentation. Despite these 
effects, and contrary to several past studies, we find that strongly-perturbed molecular cloud 
cores are able to fragment to form wide binary systems even in the presence of quite strong 
magnetic fields. For massive circumstellar discs, we find that slowing of the collapse caused 
by the magnetic support decreases the mass infall rate on to the disc and, thus, weakens gravi- 
tational instabilities in young massive circumstellar discs. This not only reduces the likelihood 
that they will fragment, but also decreases the importance of spiral density waves in providing 
angular momentum transport and in promoting planet formation. 

Key words: (magnetohydrodynamics) MHD - magnetic fields - star formation - binary stars 
- circumstellar discs 



1 INTRODUCTION 

The ability of a magnetic field to support a cloud against gravity 
may be quantified using the mass-to-flux ratio, which for a spherical 
cloud is given by 

M M 



(1) 



where M is the mass contained within the cloud volume, $ is 
the magnetic flux threading the cloud surface at radius R assum- 
ing a uniform magnetic field Bq. By a straightforward applica- 
tion of the virial theorem it can be shown that there exists a crit- 
ical value of A//$ below which a cloud will be supported against 
gravitational collapse by the magnetic field. For a u niform, spher- 
ical cloud this critical va lue takes the form (e.g. iMestell 1 19991 : 
iMac Low & Klesserj|2004l) . 
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where G and are the gravitational constant and the permeability 
of free s pace respectively and ci i s a p arameter determined numer- 
ically bv lMouschovias & SpitzeJ \\91^ to be ci ~ 0.53. 



Obse rvations of magnetic field s in molecular clouds jCrutcheJ 

1 19991 : see lHeiles & Crutchej|200^ for a recent review) show that 
magnetic fields are routinely observed with strengths where mag- 
netic and gravitational forces are comparable (i.e. the mass-to-flux 
ratio is close to the critical value). Thus, magnetic effects must be 
included in a complete theory of star formation. Despite this fact, 
most three-dimensional numerical simulations of star formation 
have neglected the effects of magnetic fields. To date, only a few 
works have begun to examine their impact on the star-formation 
process. 

The earliest three-dimensional magnetohydrodynamic (MHD) 
simulations of molecular cloud cores were those perfor med by 
iDorfil ( Il982i) using a Cartesian grid code and iBenj ( Il984l) using a 
smoothed particle hydrodynamics (SPH) code. Both studies found 
that magnetic fields inhibit the collapse of molecular cloud cores 
and could lead to the formation of bar-shaped structures during col- 
lapse. Dorfi also noted that magnetic braking can remove the bulk 
of a core's angular momentum over a period of time that is short 
compare with the core's free-fall time. Indeed, magnetic fields are 
often invoked to solve the 'angular momentum problem' in star for- 
mation by braking the rotation of molecular cloud cores. 
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IPhiUipj i ll986allbh aimed to investigate whether magnetic 
fields enhanced or inhibited fragmentation of collapsing molecu- 
lar cloud cores using an early MHD SPH code. He found that none 
of his magnetised clouds fragmented. However, given that none of 
his unmagnetised clouds fragmented either, it is likely his results 
suffered from a lack of spatial resolution due to the use of a large, 
spatially-uniform smoothing length. 

More recently. Boss performed hydrodynamical simulations 
of the collapse of molecular cloud cores where the effect of mag- 
netic fields was approximated by increasing the thermal pressure to 
mimic magnetic pressure and modelling magnetic tension forces as 
a dilution of self-gravity (Boss 2000, 2002). He found that approxi- 
mating magnetic tension effects as a dilution of self-gravity lead to 
fragmentation being enhanced. However, these calculations did not 
model the full MHD equations which may impact the applicability 
of his conclusi ons. 

Hosking & WhitworthI ( l2004ah returned to the problem of the 
impact of magnetic fields on fragmentation using an SPH code with 
variable smoothing lengths (spatial resolution). They used similar 
MHD equations to those of Phillips, but in addition they used a 
two-fluid (ion/neutral) approach that allowed them to model par- 
tially ionised fluids and take into account the effects of ambipo- 
lar diffusion. They began with a rotating molecular cloud core that 
collapsed to form a central object surrounded by a disc which frag- 
mented to form companion objects when the simulation was per- 
formed without magnetic fields. They then investigated the effect 
of magnetic fields. Beginning with a field strong enough to make 
the initial conditions subcritical, the cloud evolved via ambipo- 
lar diffusion until it became supercritical and collapsed. However, 
they found that angular momentum was transported outward due 
to magnetic braking at such a rate that fragmentation did not oc- 
cur. Indeed, in most of their cases, a disc did not even form around 
the central object; the coupling between the ions and neutrals had 
to be reduced by a factor of 25 even to obtain a small disc. Thus, 
they concluded that magnetic fields inhibit fragmentation. How- 
ever, they were also unable to follow the calculations for long peri- 
ods of time because the divergence of the magnetic field was not 
constrained to be zero and the calcul ation had to be terminated 
when it grew too large. Ziegler (200^) and ^Froman g et aL I j200d) 
each performed similar calculations to Hosking & Whitworth, but 
under the ideal MHD approximation, to test new versions of the 
NIRVANA finite-difference code and the RAMSES Godunov adap- 
tive mesh refinement code, respectively. They obtained very simi- 
lar results to those of Hosking & Whitworth both with and without 
magnetic fields (i.e., fragmentation without magnetic fields and the 
for mation of a single object when magnetic fields were included). 

iMachida et al] ( l2004l2005al lbl) followed the collapse of a wide 
range of initially magnetised rotating molecular cloud cores using 
a three-dimensional grid-based MHD code. They found that mag- 
netic fields suppressed fragmentation, but presenting their results 
in a rotation versus magnetic field plot they found that fragmen- 
tation still occurred if the cloud rotation rate was increased along 
with the magnetic field strength (i.e., for a given field strength, there 
was a rotation rate below which the cloud did not fragment and for 
stronger fields this critical rotation rate increased). Fragmentation 
in their calculations occurred via either ring fragmentation (weaker 
field strengths) or bar fragmentation (stronger field strengths). 

Finally, three-dimensional MHD simulations have also 
begun to be employed to s tu dy protostellar jet fo r mation 
jMatsumoto & T omisakal |2004| : iMachida et all l2005bl l2006l : 
iBaneriee & Pud ritz 2006). In particular, Banerjee & Pudritz per- 
formed calculations similar to those we present here. They began 



with a rotating dense core embedded in a static hot low-density 
medium with a weak magnetic field (ratio of thermal pressure to 
magnetic pressure of 84) parallel to the rotation axis of the core. 
They found that the rotating core lost angular momentum due 
to the launching of torsional Alfven waves into the low-density 
medium. The simulations began with a rotation rate that, in the 
absence of a magnetic field, resulted in the formation a ring that 
fragmented into a binary. With a magnetic field, the size of the 
ring was reduced dramatically (by approximately two orders of 
magnitude) due to angular momentum transport. As the simulation 
progressed, they observed the outward propagation of a large 
magnetically-driven bubble above and below the midplane and, 
later, the formation of a sub-AU scale jet from the central object. 

In summary, the calculations discussed above clearly confirm 
that magnetic fields act to transport angular momentum within ro- 
tating cores and to external media. Most past calculations also in- 
dicate that magnetic fields suppress the formation of binary and 
multiple systems. However, it is still not yet clear what level of 
magnetic field is required to inhibit binary formation and the depen- 
dence of this critical field strength on the initial cloud parameters 
(e.g. density perturbations, rotation rate, and/or turbulence). For ex- 
ample, Phillips ( 1986a, b) and Hosking & Whitworth (2004a) find 
no fra gmentation in any of their c alculations that included magnetic 
fields. lBaneriee & Pudritzl j2006l) produce a ring that does manage 
to fragment into a binary system , albeit one with a very small sepa- 
ration of less than 0. 1 AU, while IMachida et alj ( l2005bl) are able to 
form binaries with large separations of ~ 100 AU. 

In this paper, we report the first results from numerical cal- 
culations of the collapse of molecular cloud cores using a recently 
developed method for Smoothe d Particle Magnetohydrodynamics 
(SPMHD) dPrice & MonaghaiJ 12005, hereafter PM05). The calcu- 
lations were performed in the ideal MHD approximation. Thus, 
we are limited to studying the collapse of clouds that are ini- 
tially supercritical (i.e. they collapse without the need for ambipo- 
lar diffusion). It is straightforward to extend the method in the fu- 
ture to model ambipolar d iffusion (e.g. by adapting the two-fluid 
method of Hosking & Whitworthll2004bl) . However, given the dif- 
fering conclusions that currently exist in the literature as to the ef- 
fects of magnetic fields on fragmentation as discussed above, we 
have chosen to begin with the simplest case first. Our aims are first, 
to demonstrate the potential of this new method for star formation 
applications, and second, to investigate the issue of how magnetic 
fields impact the evolution and fragmentation of collapsing molec- 
ular cloud cores due to the extra support and the angular momentum 
transport they provide. We note in passing that a parallel but com- 
plementary approach to developing an algorithm for MHD in SPH 
has been un dertaken by Steinar B0rve and colleagues ( lB0rve et al.l 
I2OOIU2OO6I) based on a regularisation of the underlying particle 
distribution, although this algorithm has not yet been applied to 
star formation problems. 

The numerical method is discussed in ij2| The setup of the 
simulations is discussed in ijsjarid results are presented in i}4| These 
results are discussed in ^and summarised in ^ 



2 NUMERICAL METHOD 
2.1 Hydrodynamics 

The hydrodynamic method is based on a binary tree sm oothed par- 
ticle h ydrodynamics (SPH) code originally wr itten bv iBenz et al.l 
( Il990t) . Important modifications were made bv lBatel ( ll9'95h in the 
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form of sink particles (see below) and the use of individual par- 
ticle timesteps. In developing the code for MHD calculations we 
have also made several improvements to the hydrodynamics, most 
significantly in the form of incorporating the so-call e d 'variable 
smoothing length terms' jSpringel & HernquisllllOojl : iMonaghanI 
l2002l : |Pricell2004) which ensure that both energy and entropy con- 
servation are maintained to timestepping accuracy. The smoothing 
length for each particle is spatially adapted with the density accord- 
ing to the rule 




where h is the smoothing length, m is the particle mass, p is 
the density and is a dimensionless parameter which we set to 
h = 1.2, corresponding to an average of ~ 60 neighbours for 
each particle. In the variab le smoothing length formulation of 
ISpringel & HemquistI ( |2003) this is a non-linear equation for both 
h and p (that is, pi for particle i is calculated via an SPH summation 
which depends on hi) that we solve iteratively jPrice & Monaghar] 
[2OO6) . The SPH equations are integrated using a standard second- 
order leapfrog m ethod with individual timesteps for each particle 
teate etaljl99^ . 

Gravitational force softening is performed using the usual 
SPH cubic spline kernel with an adaptive softening length set equal 
to the loc al smoothing length. We u se the formalism developed re- 
cently bv lPrice & MonaghanI l l2006h which ensures that momentum 
and energy conservation are retained even in the presence of a spa- 
tially variable softening length. 

Artificial viscosity and, for MHD, resistivity are applied as in 
PM05 w hich in the case of viscosity, is equivalent to the formalism 
given bv lMonaghanI ( Il997h with the dimensionless parameters con- 
trolling the artificial viscosity and resistivity allowed to evolve with 
time as in Morris & Monaghan ( 1997) (genera lised for resistiv ity in 
PM05). It should be noted that the lMorris & MonaghanI j 19971) arti- 
ficial viscosity switch is not particularly useful during the collapse 
phase as it simply responds to the velocity divergence due to the 
gravitational infall (rather than shocks). Whilst alternative switches 
could be investigated, here we simply cap the maximum value for 
the viscosity parameter to 1.1. 



2.2 Equation of state 

We use a barotropic equation of state for the thermal pressure of 
the form 

P = Kp''. (4) 

The polytropic exponent 7 changes according to 

7 = 1, P < 10 g cm , 
7 = 7/5, p > lO^^^'g cm~^, 

giving an isothermal equation of state for low densities but a tran- 
sition at higher densities which reflects the heating of the gas as it 
becomes optically thick. The value of K is set equal to the square 
of the isothermal sound speed Cs at low densities and changes at the 
transition density so as to make the pressure a continuous function 
of density. We set the transition of the equation of state at a compar- 
atively low density (10^ ^*g cm^"^) in order to inhibit disc fragmen- 
tation in the calculations, since in this paper we are interested in the 
effects of magnetic fields on fragmentation that is seeded by initial 
density perturbations rather than on whether the discs that form in 



the system subsequently fragment to produce a higher-order multi- 
ple system. 

2.3 Sink particles 

Sink particles were introduced bv lBate et al.l ( Il995h in order to fol- 
low star formation calculations beyond the initial collapse phase in 
an efficient manner. When pressure-supported fragments form in 
the calculations due to the parametrized equation of state discussed 
above, the high densities and short dynamical times means that it 
becomes computationally expensive to follow the internal evolution 
of these protostars over timescales comparable to the initial free- 
fall time of the cloud. Instead, for the calculations presented here, a 
sink particle is inserted once the peak density exceeds ps ~ lO^^^g 
cm~^. 

The sink particle is formed by replacing the SPH gas parti- 
cles contained within race = 6.7AU of the densest gas particle in 
a pressure-supported fragment by a point mass with the same mass 
and momentum. Gas which later falls onto this point mass is ac- 
creted if it is bound and has a specific angular momentum which is 
lower than that required to form a circular orbit at the radius race 
from the sink particle. Sink particles interact with the gas only via 
gravity and accretion. When magnetic fields are present, this means 
that the magnetic field in the central regions is removed by the sink 
particle and consequently not allowed to feed back on the surround- 
ing cloud. Whilst this is obviously a somewhat crude approxima- 
tion at present, it enables us to study the subsequent accretion of 
gas to form a circumstellar disc and, for calculations that form bi- 
naries, the subsequent evolution of the binary system without the 
calculation grinding to a halt trying to follow the internal dynamics 
of the protostars themselves. 

2.4 Magnetohydrodynamics 

The magnetohydrodynamics implemented is an extension 
of the methods fo r MH D in SPH developed recently by 
IPrice & Monaghan' llOOAa'h) and PM05. The extension is 
that here, as in Price & R osswog ( 2006) , the divergence constraint 
on the magnetic field is satisfied by expressing the magnetic field 
as a function of two scalar variables a and /3 according to 

B = Va X V/3, (5) 

which we refer to as the 'Euler potentials ' jStemI I1970I) 
but elsewhere referred t o as the 'Clebsch formalism' 
jPhilUps & Monaghad Il985h or 'flux co-ordinates'. We cal- 
culate the gradients of the potent ials using a n SPH summation 
which is exact for a linear gradient ( |Pricel2004h . that is, 

Xmi'V"! = - ^ - aj)V^IVij(/ii), (6) 

J 

= -^m,(ft-/3,)Vrmi(/iO, (7) 
J 

where latin indices refer to the particles, greek indices refer to vec- 
tor components and the matrix quantity Xmi^ is given by 

XM. = E"^j(^''-'^i)^"'^»:'('^0- (8) 
j 

The gradient calculations involve solving a 3 x 3 matrix equation 
in each case, but this is straightforward to solve analytically and 
implementation in the code requires only the temporary storage of 
the matrix variable x for a given particle (ie. we do not need to 
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store calculated values for more than one particle). Furthermore, 
since the summations l|7j and ^ do not involve the density on 
neighbouring particles, the gradients can be calculated efficiently 
alongside the usual density summation. The B field calculated from 
the Euler potential gradients in this manner is then used in the force 
equation as in PM05. 

Apart from satisfying the divergence constraint, the other main 
advantage of using the Euler potentials representation is that, for 
ideal MHD, they evolve according to 



(9) 



corresponding to t he advection of magnetic field lines by La- 
grangian particles l lStern|[l969) . We extend the Euler potentials 
method to non-ideal MHD by incorporating shock-capturing dis- 
sipation terms, the form of vy hich is given, using a simple generali- 
sation of the terms derived in lPrice & Monagharj j2004bl) . by 



= Y^m,^^^{a,~a,)\VWrj\, (10) 

diss j 

^ Y^rn,^^{P.~l5,)\VW-\, (11) 

diss j P^j 



where the summation is over neighbouring particles, Wsig is a maxi- 
mum signal velocity between the particle pair as in PM05, the mean 
density p = 0.5(pi -I- pj), \S/Wij\ refers to the magnitude of the 
mean kernel gradient \/Wij — 0.5[\/Wij{hi) + \7Wij{hj)] and 
is a time-variable co-efficient for each particle that is evolved 

as described in PM05. 

The magnetic force is computed using the lMorrisl ( ll996iJ) for- 
malism discussed in PM05 (using the B computed from equation[5]l 
which ensures stability of the SPMHD formalism against particle- 
clumping instabilities in the regime where gas pressure is dominant 
over magnetic pressure. The MHD part of the force equation (that 
is, apart from the gravitational and artificial viscosity forces) reads 



dt 



Pi + ^Bf 
^•■tPi 

{B^B"), - {B''B'')j 



Figure 1. Initial setup. A uniform density cloud of radius 4 X lO^^cm and 
mass 1 Mq in uniform rotation is placed in a warm, low density surrounding 
medium (density ratio 30:1). Magnetic field lines thread the cloud. For the 
binary calculations an m = 2 perturbation is applied to the cloud density. 
The figure shows a cross section slice of the density at j; = together with 
magnetic field lines aligned with the rotation (z) axis. 



3 INITIAL CONDITIONS 

The initial cloud is a sphere of radius i? = 4 x 10^^ cm (0.013 
pc) and mass M — IMq with mean density po — 7.43 x 
lO^^^g cm"^. The free-fall time of the cloud is tg = 2.4 x 10'' 
years. 

We assume, for simplicity, that an initially uniform magnetic 
flux threads the cloud and connects it to the surrounding interstel- 
lar medium. However, a key factor in the problems studied here is 
the angular momentum transfer introduced by the magnetic field 
in the form of magnetic braking of the rotating core. Thus careful 
attention must be paid to the boundary condition at r = i?. Ex- 
periments with simple boundary conditions for SPH (for example, 
using constant pressure boundaries or ghost particles) proved some- 
what unsatisfactory, particularly because, in the higher magnetic 
field strength runs, significant material is flung outwards by the 
cloud along the magnetic field lines into the surrounding medium. 
We therefore model the boundaries self-consistently by placing the 
cloud within a uniform, low density box of surro unding material 
in pressure equilibrium w ith the cloud (see also iHoskindliooi: 
Baneriee & Pudritzl|200 i). To ensure regularity of the particle dis- 



poptpj 



where Vi , pi , Pi and Bi refer to the density, pressure and mag- 
netic field of particle i and f2 is a normalisati on term related to 
the gradient of the smoot hing length (as in, e.g. lMonaghar]|2002l : 
IPrice & Monagharj 1200^. The first term in l ll2b is the isotropic 
hydrodynamic + magnetic pressure force and the second term is 
the magnetic tension force. For simulations "without magnetic ten- 
sion" we do not include the latter term. 

A detailed summary of the recent changes to the hydrody- 
namic method (including details of the implementation of the vari- 
able smoothing length SPH formalisms in both the pressure and 
gravity terms) and the implementation of the Euler potentials into 
the numerical code (including test problems comparing the use of 
them to th e 'standard' SPMHD form alism of PM05) are discussed 
in detail in lRosswog & Pric3 ( |2007|) and we refer the reader to this 
paper for an up-to-date summary of the present numerical code (the 
specific code described differs from that used here but the algo- 
rithms implemented in each are identical). 



Q2y"ibutions at the box boundaries, we use quasi-periodic boundary 
conditions at the box edge (that is, particles within 2 smoothing 
lengths of the boundary are 'ghosted' to the opposite boundary, 
with no self gravity between SPH particles and ghost particles). 
Since a uniform magnetic field necessarily implies a linear gradi- 
ent in the Euler potentials, continuity of the magnetic field across 
the box boundary is ensured by adding an offset to the values of a 
and P copied to the ghosted particles corresponding to an extrapo- 
lation of the linear gradients outside the box boundaries. 

We find that satisfactory results are obtained using a box size 
of —8 X 10^*'cm < x,y, z < 8 X 10^*'cm (that is, twice the cloud 
radius in each direction) and a density ratio of 30 : 1 between the 
cloud and the surrounding medium. This density ratio was chosen 
simply to ensure that the surrounding medium is sufficiently hot so 
as not to contribute significantly to the self-gravity of the cloud (that 
is, Cg > GAI/R). The initial setup is shown in Figure[T] showing a 
cross-section slice of density aty — with overlaid magnetic field 
lines for a field initially oriented in the direction. 

Both the spherical cloud and the surrounding medium are set 
up by placing the particles in a regular close-packed lattice arrange- 
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ment (e.g. Hoskind 2002') which is a stable arrangement for the 
particles jMorrij 1996b). Whilst such regularity introduces some 
undesirable side effects due to the lattice regularity in the initial 
collapse phase, these small and transient effects are largely elimi- 
nated by the time star formation occurs. The exact positions of the 
box boundaries are adjusted slightly to ensure continuity of the lat- 
tice across the ghosted boundary. 

The magnetic field strength is characterised in terms of a 
mass-to-flux ratio expressed in units of the critical value l|2](. The 
corresponding magnetic field strength is given by 



Bo = 814^G 



M 



R 



O.OlSpc 



(13) 



where M/"I> is the mass to flux ratio in units of the critical value. 

For both the axisymmetric and binary star formation prob- 
lems, we calculate a sequence of collapse models with mag- 
netic field strengths corresponding to mass-to-flux ratios, in 
units of the critical value (O, of oo (that is, hydrodynamic), 
100,20,10,7.5,5,4,3,2 and 1. For our choice of cloud mass 
and radius these correspond to field strengths of B — 
0,8.1,40.7,81.3,108.5,163,203,271,407 and 814^iG' respec- 
tively. The initial Alfven speeds in the cloud are therefore 0, 8.4 x 
10"^0.42, 0.84, 1.12, 1.68, 2.10, 2.80,4.21 and 8.42 km/s re- 
spectively. We compute these sequences in each case both for a 
magnetic field initially aligned with the rotation axis (ie. oriented 
in the direction) and for a magnetic field initially oriented in the 
a:— direction, thus bracketing two geometric extremes. We find that, 
whilst there are similar global trends in both cases due to the effect 
of magnetic pressure, the orientation of the magnetic field with re- 
spect to the rotation axis plays an important role in determining the 
outcome. 

Simulations are performed using 300,000 particles in the 
cloud itself. Including the external medium as discussed above, this 
results in a total of 451,233 particles in each simulation. We have 
also performed simulations using 30,000 and 100,000 particles in 
the cloud, which show only minor differences in results (see dis- 
cussion in ij4.2| l. In order to resolve the local Jeans mass throughout 
the calculations, we require at least 30 ,000 particles for our chosen 
equation of state jBate & Burkert|[l997i) . Thus, all our calculations 
resolve the local Jean mass, and the high resolution calculations 
presented here do so by an order of magnitude in particle number. 



4 RESULTS 

4.1 Axisymmetric collapse 

For an axisymmetric collapse, we set the initial cloud to be uniform 
density in uniform rotation with an angular velocity ofQ — 1.77 x 
10"^"^ rad s~^ corresponding to a ratio of rotational to gravitational 
energy /3r ~ 0.005 and Qta = 0.136. 

The initial cloud temperature is set such that the ratio of ther- 
mal to gravitational energy a = 0.35 (where by gravitational en- 
ergy we mean the magnitude of the gravitational potential energy). 
This corresponds to an internal energy of 7.04 x lO^ergs g~^ and, 
assuming a mean molecular weight for molecular hydrogen (ie. 
H = 2), corresponds to an isothermal temperature of 11. 3K. The 
initial sound speed in the cloud is Cs — 2.16 x 10^ cm s~^ (for 
comparison the sound speed in the external medium is Cs, medium = 
11.9 X 10* cm s-\ ie. c^ ^cdium - 4.2GA//i?). 

Expressed in terms of the ratio of gas to magnetic pres- 
sure, /3m, our chosen sequence of mass to flux ratios of 




Figure 2. Column density and projected magnetic field lines in the collaps- 
ing, axisymmetric cloud, showing results at iff = 1.01 using mass-to-flux 
ratios in units of the critical value of (left) 100 (that is, a very weak field) 
and (right) 3 (ie., a very strong field) with an initial field aligned with the 
rotation (z) axis. In the weak field case (left) the magnetic field lines are 
deformed by the cloud, leading to an almost-spherical collapse, whereas 
for strong fields (right) inaterial is strongly channelled along magnetic field 
lines, leading to an bipolar-shaped cloud and a pancake-like collapse. 



cxD, 100, 20, 10, 7.5, 5, 4, 3, 2 and 1 in this cloud correspond to 
/3m = oo, 983, 39, 9.8, 5.5, 2.5, 1.6, 0.85, 0.39 and 0.098, respec- 
tively. 

We compute these sequences for initial magnetic fields thread- 
ing the cloud which are aligned parallel and perpendicular to the 
rotation axis (that is, in the z— and x— directions respectively in 
our computational domain). 



4.1.1 Initial field aligned with the rotation axis 

As an illustration of the global cloud evolution for runs with the 
field aligned with the rotation axis, Figure |2] shows the column 
density and projected magnetic field lines in the collapsing cloud 
at t = 1.01 free fall times for two runs with a very weak (A//<I> = 
100, left panel) and very strong (A//^ = 3, right panel) mag- 
netic field, projected in the z ~ x direction (integrated through the 
direction). For runs with weak magnetic fields (Af/$ > 10), 
the collapse is almost spherical and the magnetic field lines are 
bent by the infalling gas flow, developing a toroidal component of 
similar magnitude to the field. For very strong magnetic fields 
(Af/<I> < 3) the collapse is strongly channelled along the magnetic 
field lines, producing a pancake-like collapse and developing only 
a small toroidal component in the large scale magnetic field. In the 
latter case the interaction between the strong magnetic field and the 
rotating cloud is also evident by the shape of the cloud in the right 
hand panel of Figure (2] which has evolved to a bipolar configura- 
tion as material is flung outwards along the magnetic field lines. 

The magnetic field strength in the collapsing cloud is found 
to scale with density approximately as B oc p" '', with a pro- 
nounced flattening off in the disc itself for the higher field strength 
runs. This is in good agr eement with a similar result found by 
ISaneriee & Pudrita ( l2006h who considered only a collapse with a 
relatively low field strength, as here with the initial field aligned 
with the rotation axis. 

The results from axisymmetric collapse calculations using an 
initial field aligned with the rotation axis are shown in Figure [3] 
showing column density through the cloud (integrated through the 
2;— direction) at 5 different times (left to right) for a sequence of 
runs of increasing magnetic field strength (top to bottom), where 
time is shown in units of the free-fall time and magnetic field 
strengths are expressed in terms of the mass-to-flux ratio in units 
of the critical value. The figure shows the runs with A//<I> = 
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Figure 3. Results of the axisymmetric collapse calculations with initial magnetic field aligned with the rotation (z) axis, showing column density in the 
collapsed cloud (integrated through the z direction). Columns from left to right show snapshots at a given time (given in units of % = 2.4 X 10* yrs), whilst 
from top to bottom rows show the results for increasing magnetic field strengths, given as a mass to flux ratio in units of the critical value but which correspond 
to S = 0, 40.7, 81.3, 108.5, 163 and 203^G respectively. Increasing the magnetic field strength tends to delay and also suppress disc formation. 
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Time [ free fall times ] 

Figure 4. Circumstellai' disc mass (top) and radius (bottom) plotted as a 
function of time in tlie axisymmetric calculations using a field aligned with 
the rotation axis (see Figure|3). The disc mass is defined simply as the mass 
of material which has a density > 10~^''g/cm^ (ie. an order of magnitude 
above where the gas becomes optically thick) and the disc radius is defined 
as the radius containing 99% of the mass of the star+disc system. For veiy 
low magnetic field strengths the disc is more massive and thus dominated 
by gravitational instabilities resulting in spiral arms, producing a disc radius 
which is oscillatory in time. Increasing the magnetic field strength tends to 
suppress disc formation, resulting in smaller, less massive discs which are 
less dominated by self-gravity. 



00,20,10,7.5,5,4 and 3. The M/$ = 100 run differs only 
slightly from the hydrodynamics case and has not therefore been 
plotted. For very strong magnetic fields {M /$ < 3) the collapse 
is strongly inhibited as the mass-to-flux ratio is close to the criti- 
cal value. Gas is strongly channelled along the magnetic field lines 
(Figure |2ll and we find that a ring is formed at the cloud radius 
(ie. far away from the central regions, and long before any collapse 
has occurred in the centre) in the mid-plane which becomes grav- 
itationally unstable and fragments. Because this fragmentation is 
essentially a boundary effect, we do not discuss these calculations 
any further. 

The results for intermediate field strengths show a clear trend 
in both the formation of the protostar and the subsequent size of 
the disc which forms (Figure [3}, namely that protostar formation 
occurs progressively later as the field strength is increased and the 
disc which forms is smaller, less massive and thus also less prone 
to gravitational instabilities such as spiral arms which are clearly 
evident in the hydrodynamic run (top row). 

The delay in the onset of fragmentation is largely a result of 
the extra support provided to the cloud by the magnetic pressure. 
The size and mass of the resultant disc is affected in this case 
mainly by the subsequently lower accretion rate onto the central 
core. Magnetic tension is found to play only a relatively minor role 



when the field is aligned with the rotation axis. The role of mag- 
netic pressure versus magnetic tension is apparent in the compar- 
ison between runs using fields aligned and fields perpendicular to 
the rotation axis, discussed below. 

We have quantified the results shown in Figure [3] by calculat- 
ing the disc mass and radius as a function of time, plotted in Fig- 
ure|4] We define the disc mass as simply the mass of material which 
has a density > 10~^''g/cm'^ (ie. an order of magnitude above the 
density at which the gas has become non-isothermal according to 
our equation of state[5j- Similarly the disc radius is defined as the 
radius which contains 99% of the mass of the star-l-disc system. 
For the lower field strength runs (e.g. the solid line in Figure |4] 
corresponding to the hydrodynamic run), a disc is formed with a 
radius of ~ 30 — 50 AU at around 1 free-fall time. The disc ra- 
dius shows pronounced oscillations in time because the relatively 
massive discs are dominated by gravitational instabilities which 
drive non-axisymmetric features such as spiral arms. Increasing the 
magnetic field strength tends to suppress disc formation, producing 
smaller, less massive discs with masses and radii which increase 
more slowly in time and which show a smoother evolution due to 
the suppression of gravitational instabilities. 

We have also performed axisymmetric calculations with 
fields using higher rotation rates (/3r = 0.01 and [3^ — 0.16, not 
shown). Increasing the rotation rate tends to enhance gravitational 
instabilities in the disc and thus the ability to fragment into multi- 
ple systems, however we find that the effect of the magnetic field is 
the same - namely that increasing the magnetic field strength sup- 
presses gravitational instabilities in the disc and thus the propensity 
to fragment. 



4.1.2 Initial field perpendicular to the rotation axis 

The results of axisymmetric calculations starting from a magnetic 
field initially perpendicular to the rotation axis (that is, a field ini- 
tially in the .t— direction) are shown in Figure |5] shown at the 
same times as in Figure [3] for runs using mass-to-flux ratios of 
20, 10, 7.5, 5 and 4 (top to bottom), corresponding to the middle 
5 rows of Figure [3] As in the aligned-field calculations increas- 
ing the magnetic field strength can be seen to suppress the forma- 
tion of the circumstellar disc, forming smaller, less massive discs 
in which the spiral structure produced by gravitational instabilities 
is suppressed. However in this case the geometry of the magnetic 
field results in an increasingly bar-like collapse (most evident in the 
M/<E> = 4 and Af/<E> = 5 runs at ts = 1.06 in Figure^ which 
is wound up by the rotation of the core to produce a low density 
spiral structure surrounding the central disc. The winding up of the 
bar results in oppositely directed field in the spiral structure which 
produces significant reconnection, particularly in the stronger field 
runs where this reconnection results in a small amount of ejected 
material from the central regions. 

The evolution of disc mass and radius is similar to the aligned- 
field case (see Figure |4ll, mainly because the changes due to mag- 
netic tension are in the lower density material surrounding the disc 
rather than in the disc itself. However, we note that the misaligned 
field seems to inhibit gravitational instabilities in the dense central 
part of the disc even more than in the aligned-field case. 

4.2 Binary star formation 

The binary star formation problem we consider is a variation o n 
the 'standard isothermal test case ' of lBoss & Bodenheimeil ( Il979t) . 
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Figure 5. Results of the axisymmetric collapse calculations with initial magnetic field perpendicular to the rotation axis, which may be compared with the 
aligned-field runs shown in Figure[3] As in the aligned field case, increasing the magnetic field strength tends to delay and also suppress disc formation. In this 
case, however, the magnetic field geometry results in a increasingly bar-like collapse. For lower field strengths this bar-like structure is wound up to produce 
low density spiral structure surrounding the disc. 



the main differences in the hydrodynamic problem being the larger 
cloud radius. The initially uniform density cloud is given a non- 
axisymmetric m = 2 density perturbation with amphtude A of the 
form 



p = po[l + Acqs {24>)], 



(14) 



where (j> is the azimuthal angle about the rotation (z) axis and po 
is set as previously. The cloud is in uniform rotation with Q. — 
1.006 X 10~^^ rad s~^, corresponding to a ratio of rotational energy 
to gravitational potential energy of /3r — 0.16 and Qtg = 0.77. 
The internal energy of the cloud is set such that the ratio of thermal 
to gravitational energy is q = 0.26. This corresponds to an inter- 
nal energy of 5.23 x lO^ergs g~^, a sound speed in the cloud of 
1.87 X lO^cm/s and, assuming a mean molecular weight of p = 2, 
an isothermal temperature of 8.4K. Given the 30: 1 density ratio be- 



tween the cloud and the surrounding medium and the condition of 
pressure equilibrium between the two, this means that the sound 
speed in the external medium is Cs, medium ~ 10.2 x 10'*cm/s, and 

<mcdium/(GA//i?) = 3.1. 

The initial density perturbations with amplitudes A = 0.1, 0.2 
were applied, retaining equal particle masses, by perturbing the 
particle positions from a uniform distribution according to the lin- 
earized continuity equation 



V 



(Sr) = ^. 
Po 



(15) 



Using il4\ the perturbation in position is given by 

S<t>^ sin {2<l>o), (16) 
where (f>o is the unperturbed azimuthal angle. For calculations using 



The impact of magnetic fields on single and binary star formation 9 




Figure 6. Results of the binaiy stai' formation calculations using a magnetic field aligned with the rotation axis. As in Figure [3] column density in the 
fragmenting cloud is plotted, in a sequence of snapshots in time (from left to right) and increasing magnetic field strength (from top to bottom). As previously 
times are given in units of the free-fall time, tff = 2.4 X 10* yrs and magnetic field strength is expressed in terms of the mass-to-flux ratio in units of the 
critical value, corresponding to B = 0, 40.7, 81.3, 108.5, 163 and 203/iG from top to bottom respectively. The separation of the binary system decreases 
with increasing magnetic field strength, eventually forming only a single star sun'ounded by a rotationally-dominated circumstellar disc. 



A > 0.5 this is no longer a good approximation an we instead apply 4.2.1 Initial field aligned with the rotation axis 
the perturbation by changing particle masses appropriately. 

The results of the binary star formation calculations with the initial 
magnetic field aligned with the rotation axis are shown in Figure|6] 
where as previously we have computed a series of runs of increas- 
ing magnetic field strength, corresponding to mass-to-flux ratios in 
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units of the critical value of (from top to bottom) oo (that is, hydro- 
dynamics), 20, 10, 7.5, 5 and 4. Given the initial temperature and 
density of the cloud, the coiTesponding values for the ratio of gas 
to magnetic pressure are /Jm = oo, 39, 9.8, 5.5, 2.4 and 1.57 re- 
spectively. The results are shown at 5 different times (left to right) 
given in units of the free-fall time (2.4 x lO^yrs). The global evolu- 
tion of the cloud prior to star formation is similar to that discussed 
in the axisymmetric case (see i]4.1l above). Collapse in the runs 
with M/<1> < 3 is strongly suppressed by the magnetic field and 
these runs are therefore not shown. 

In the binary case there is a clear trend of delayed collapse 
and decreasing binary separation as the magnetic field strength in- 
creases (top to bottom in Figure|6}. Eventually, for very strong mag- 
netic fields (7l//$ < 4) binary formation is inhibited completely 
(the remnant of the initial m = 2 density perturbation is just visi- 
ble at tff — 1.3 in the M/$ = 4 case, but merges to form a single 
protostar at iff = 1.33 and is subsequently surrounded by a mas- 
sive, circumstellar disc). We follow the calculations to the point at 
which sub-fragmentation occurs in the weak magnetic field runs 
(ie. at tfi = 1.46 for the hydrodynamic and M/$ = 20 runs). 

Whilst it is tempting to attribute the decrease in binary separa- 
tions observed in Figure[6]to angular momentum transport induced 
by the magnetic field (ie. via magnetic braking of the core), this is 
in fact not the dominant cause. Rather, the magnetic pressure acts to 
increase the effective "a" (ratio of thermal to gravitational energy) 
in the cloud which reduces the propensity of the cloud to fragment. 
That this is the case is demonstrated in Figure|7]where we show the 
results of three runs with M/^ = 7.5,5 and 4 that were performed 
a) without magnetic tension forces (top row of Figure [7J and b) 
with no magnetic fields, but with the cloud temperature increased 
to give equivalent effective values of a (where the corresponding 
values are a = 0.29, 0.34 and 0.38 respectively, determined by 
transferring the initial magnetic energy to the thermal energy of the 
cloud). The results are shown at ts = 1.33, corresponding to the 
centre panels of the fourth, fifth and sixth rows in Figure|6] The in- 
teresting point to note is that the trend in binary separation is similar 
in both cases (there is a transition from a binary to a single proto- 
star), indicating that it is the magnetic pressure which is providing 
the dominant effect in suppressing binary formation and that the 
effect of magnetic pressure is similar to an increase in the effective 
thermal energy of the cloud. That is not to say that magnetic tension 
forces are without effect as there are differences between the results 
with and without magnetic tension forces. Nor is the effect of mag- 
netic pressure exactly equivalent to a thermal pressure, indicated 
by the slight differences which remain between the hydrodynam- 
ical calculations and the magnetic pressure-only calculations and 
the dependence on field geometry (see i]4.2.2l below). However, 
the dominant effect observed in Figure[6]is primarily attributable to 
magnetic pressure effects rather than magnetic braking. 

The results shown in Figure|6]are quantified in Figure[8]which 
shows the binary separation as a function of time for the magnetic 
field strengths shown in Figure [6] Prior to sink particle formation, 
we define the binary separation as the distance between the two 
highest density maxima constrained to be in opposite hemispheres. 
Some spurious effects from this definition are visible at tff ^ 1.1 
in the weak magnetic field runs, where there is a transient density 
maximum due to particle gridding effects. Also the results become 
meaningless when only one star is formed, as in the M/<^ = 4 case 
beyond ta ~ 1.3 (lowest solid line) at which point we do not plot 
any more points. However, the general trend is clear - namely that 
the binary separation decreases monotonically as the magnetic field 
strength is increased. The M/$ = 5 case can be seen to form a 




Figure 7. Results of (top row) simulations using a magnetic field initially 
aligned with the rotation axis but without magnetic tension forces and (bot- 
tom) purely hydrodynamical simulations using ratios of thermal to gravi- 
tational energy, a equivalent to the effective support provided by both the 
thermal and magnetic pressures. The plots are shown at tff = 1.33 for sim- 
ulations with mass to flux ratios (top row) of (from left to right) 7.5, 5 and 
4, corresponding to the centre panel of the last three rows of Figure |6] In the 
hydrodynamical case (bottom row) we have used a = 0.29, 0.34 and 0.38 
respectively. In both cases the same transition from a binary to a single pro- 
tostar is observed, indicating that, for fields aligned with the rotation axis, 
magnetic pressure plays the dominant role in suppressing binary formation 
and that this is similar to an increase in the effective thermal energy of the 
cores. 
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Figure 8. Separation of the binary systems shown in Figure |6] plotted as a 
function of time, with field strengths (mass to flux ratios) indicated by the 
legend. Prior to sink particle formation, we define the binary separation as 
the distance between the two highest density maxima constrained to be in 
opposite hemispheres. The plot quantifies the decrease in binary separation 
with increasing magnetic field strength seen in Figure |6] The dip at % 
1.1 in the low magnetic field strength runs is a spurious feature due to a 
transient density peak related to particle gridding effects. 



very close but eccentric binary system which goes through a closest 
approach of ~ 30 AU. 

For the binary star formation case we have also performed 
simulations which start from larger initial perturbations (A — 0.2 
and A — 0.5 percent) and with larger and smaller cloud radii 
(5 X 10 cm and3xl0"cm). The runs with larger initial perturba- 
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Figure 9. Results of binary star formation calculations using a magnetic field initially oriented perpendicular to the rotation axis (ie. initial field in the 
a;— direction). As previously times are given in units of the free-fall time, = 2.4 X 10* yrs and magnetic field strength is expressed in terms of the 
mass-to-flux ratio in units of the critical value, corresponding to B = 40.7, 81.3, 108.5, 163 and 203/.tG from top to bottom respectively. The transition from 
a binary to a single protostar occurs at a lower field strength than with the initial field aligned with the rotation axis. 



tions are progressively less influenced by the magnetic field. Indeed 
for the A = 0.5 percent perturbation all of the calculations up to 
a mass to flux ratio of Af/"I> ~ 2 formed binaries with separations 
that are not significantly influenced by the magnetic field strength. 
The effect of changing the cloud radius is to change at what stage 
during the coflapse the cloud becomes optically thick (ie. the equa- 
tion of state changes from isothermal to non-isothermal). Thus disc 
fragmentation is (relative to the calculations presented here) sup- 
pressed for smaller cloud radii and enhanced for larger radii but 
with similar trends in the influence of the magnetic field. 

4.2.2 Initial field perpendicular to the rotation axis 

Results of binary star formation calculations beginning with a mag- 
netic field oriented perpendicular to the rotation axis (that is, with a 



field initially in the x— direction) are shown in Figure[9] As in Fig- 
ure |6] some general trends are clear: increasing the magnetic field 
strength leads to a delayed collapse and increasingly suppresses bi- 
nary formation. In this case, however, the transition from a binary 
to a single star occurs earlier (that is, a single star is formed at 
= 5 in Figure [9] compared to Af/cl? = 4 in Figure |6} and 
the binary perturbation is increasingly deformed by the magnetic 
field, which at higher field strengths results in a "double bar-like 
collapse" (most evident in the higher field strength runs in Figure[9] 

As previously, the global trends (delayed collapse and tran- 
sition to a single star) are the result of the extra support provided 
to the cloud by magnetic pressure alone. This is demonstrated by 
Figure [To] which shows the results of similar calculations (that is, 
with fields initially perpendicular to the rotation axis) but with mag- 
netic tension forces turned off. In this case the transition to a single 
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Figure 10. Results of binary star formation calculations using a magnetic field initially oriented perpendicular to the rotation axis (ie. initial field in the 
a;— direction), as in Figure|9]but with magnetic tension forces turned off. The transition from a binary to a single protostar occurs at a lower field strength when 
magnetic tension forces are excluded (i.e., magnetic tension aids binary formation). 



star occurs for even lower magnetic field strengths (at M/$ = 10). 
This indicates not only that magnetic pressure is providing the dom- 
inant role in suppressing fragmentation but also that magnetic ten- 
sion can act to dilute the effect of magnetic pressure, even aiding 
binary formation. We note that Boss (2000, 2002) similarly con- 
cluded that magnetic tension forces can act to promote fragmenta- 
tion, albeit using a simplistic approximation to model the affect of 
magnetic fields. 

Comparing the middle column of Figure [TO] to the pressure- 
only calculations with an aligned magnetic field shown in Figure |7] 
also demonstrates that the effect of magnetic pressure is dependent 
on the field geometry, acting more like a equivalent thermal pres- 
sure when the field is aligned with the rotation axis. 

The deformation of the binary perturbation in Figure|9]is not 
evident in the tension-free calculations (Figure [TOt, indicating (as 
one might expect) that this effect is the result of the gas being 
squeezed by the magnetic field lines. It is this squeezing due to 
magnetic tension that acts to hold up the rapid transition to a single 
star with increasing field strength observed in the tension-free runs 
and thus dilute the effect of magnetic pressure in suppressing frag- 
mentation. The magnetic field, being in this case aligned along the 
binary perturbation, effectively acts as a "cushion" between the two 
stars which prevents their merging. This effect, which we hence- 
forth refer to as "magnetic cushioning", is graphically illustrated 
in Figure [TT] which shows the magnetic field (arrows in left panel, 
overlaid on a column density plot) and integrated magnetic pressure 
(right panel) in the M/<^ = 10 run (corresponding to the second 
row of Figure|9]l at ts = 1.35. The "cushion" formed by the mag- 
netic field between the two stars is clearly evident, and it is this 
"magnetic cushion" which prevents the binary system from merg- 
ing to form a single star (and also produces the wonderful symme- 
try in the spiral arms). 



The results shown in Figure |9] are quantified in Figure [12] 
which shows the binary separation as a function of time for the 
magnetic field strengths shown in Figure |9] and may be compared 
with the corresponding figure (Figure [Sjl for the runs with the field 
aligned with the rotation axis. As previously, prior to sink particle 
formation, we compute the separation of two density maxima in 
opposite hemispheres. In the stronger field runs (M/$ = 4 and 
M/$ = 5) the binary perturbation is strongly deformed by the 
magnetic field, producing the observed increase in separation ob- 
served at tfi ~ 1.2. The binary separations are in each case smaller 
than the equivalent runs using a field aligned with the rotation axis, 
which demonstrates that the effect of the magnetic field of the bi- 
nary system is stronger in this case. The effect of magnetic cushion- 
ing is also apparent in the fact that the runs with i\f /<!> — 20, 10 and 
7.5 show a trend of increasing binary separation at closest approach 
(iff 1.35), in contrast to Figure [8] (although all the separations 
are smaller than in the aligned-field runs). 



5 DISCUSSION 

We have conducted a study of how magnetic fields affect the col- 
lapse of homogenous molecular cloud cores and cores with initial 
m = 2 density perturbations. In both cases, the presence of a mag- 
netic field produces a delayed collapsed, with a longer delay for 
stronger fields. This affect is easily attributed to the effect of the 
magnetic pressure on the collapse. The magnetic field gives extra 
support to the cloud over the thermal pressure alone; rather than act- 
ing like a cloud whose ratio of thermal energy to the magnitude of 
gravitational energy a — 0.26 (or 0.35 in the axisymmetric mod- 
els), the effect of the magnetic field is to raise the effective value of 
a. 
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Figure 11. Magnetic cushioning in action. The panels show column density and integrated magnetic field vectors (left) and integrated magnetic pressure 
(right) at iff = 1.35 in the Af/<I> = 10 run (coiresponding to the second row of Figure|9). The magnetic field, initially in the orbital plane, is wound up by 
the differentially rotating cloud to form a "cushion" between the binary, preventing it merging into a single protostar Thus, the magnetic cushion aids binary 
formation. 
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Figure 12. Separation of the binary systems shown in Figure |9]plotted as a 
function of time, with field strengths (mass to flux ratios) indicated by the 
legend. 



5.1 The effect of magnetic fields on protostellar discs 

In the homogenous simulations, we find that a single protostar 
(sink particle) is formed and surrounded by a disc. Stronger mag- 
netic fields lead to a delay in the formation time of the protostar 
as mentioned above, but they also decrease the rate of accretion 
onto the disc. The disc radius also increases more slowly with 
time. It is well known that the rate of infall of mass onto a massive 
disc is crucial in generating gravitational instabilities (e.g. Bonnell 
1 1994 IWhitworth et"ai]| 19951 ; iHennebelle et al.ll2004l) and, indeed, 
we see this effect here. In the purely hydrodynamical case, the disc 
surrounding the protostar is gravitationally unstable and exhibits 
strong spiral density waves soon after the protostar forms (although 



the instability is not strong enough to force the disc to fragment). 
With a magnetic field initially aligned with the rotation axis, the 
slower rate of mass infall onto the disc leads to a weakening of the 
gravitational instability such that for mass to flux ratios less than 
A//$ ~ 10 the spiral features are very weak (Figure [3l(. With a 
field initially perpendicular to the rotation axis, the gravitational 
instability is very weak even for M/^ = 20 (Figure|5](. 

Gravitational instabilities in protostellar discs may be impor- 
tant for several reasons. First, if the gravitational instability is 
strong enough , the disc may fragment to form a companion (e.g . 
Bonnelll 1 1994 i Bonnell & Bate! Il994allbl IWhitworth et"d] 1199^ : 
Rice et al. I l2005lr This is particularly rele vant to the magnetised 
star for mation simulations performed by iHosking & WhitworthI 
( l2004ij) . They began with a rotating cloud that, in the absence of 
magnetic fields formed a single object surrounded by a gravitation- 
ally unstable disc that fragmented to form companions. With mag- 
netic fields initially aligned with the rotation axis, they found that 
the disc was much smaller and did not fragment. This is consistent 
with our simulations in that we also find that magnetic fields re- 
duce the tendency for a disc to be gravitational unstable. However, 
Hosking & Whitworth attributed the inhibiting of fragmentation to 
the loss of angular momentum due to magnetic tension forces. Here 
we find that the effect of magnetic pressure in decreasing the mass 
infall rate on to the disc may be just as important in suppressing 
disc fragmentation. 

Second, if a protostellar disc is gravitationally unstable (but 
not so strongly as to fragment), spiral density waves are likely to be 
the dominant angular momentum transport mechanism within the 
accretion disc (Lodato & Rice 2004, 2005; Fromang et al. 2004). It 
is only later when the disc becomes more stable that the magne- 
torotational instability is likely to take over as the dominant angu- 
lar momentum transport mechanism. Our findings suggest that for 
magnetised cores the phase of rapid angular momentum transport 
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due to gravitational instabilities will be less important than for stars 
formed in unmagnetised cores. 

Third, it has been suggested that the spiral density waves 
in gravitationally unstable protoplanetary discs might aid the for- 
mation of large (kilometre-sized) planetesimals, and thus planets, 
by collecting centimetre and metre-sized planetesimals in the spi- 
ral waves until the density of solids along the arms is sufficient 
for gravitational instability to form solid kilometre-sized objects 
jRice et aljioosh . However, the effect that we observe here of mag- 
netic fields slowing the infall of gas onto the disc and therefore 
weakening the gravitational instability argues that these processes 
may be less important for stars forming in magnetised molecular 
cloud cores. On the other hand, it is plausible that the spiral struc- 
ture generated by the magnetic field in the outer parts of the discs 
formed in the misaligned field calculations may act in a similar 
manner 

5.2 The effect of magnetic fields on fragmentation 

Several past studies have cast doubt on the ability of strongly mag- 
netised mo lecular cloud c or es to fo rm anything other than sin- 
gle objects ( |Phillips|[T9863lhl;lHosking & Whitworth 2004a) . Even 
MHD calculations that have obtained binary systems have either 
produced binaries with very small separations jSaneriee & Pudrit3 
[2006) or have requir ed strongly magnetised cores to have very large 
initial rotation rates jMachida et ^120051^ . Given that binary and 
multiple stellar systems are very common, this is potentially a seri- 
ous problem. 

Perhaps the easiest way to obtain a binary system is to be- 
gin with a strongly perturbed initial core - a picture that is also in 
qualitative agreement with observations ( Pringlelll989h . Therefore, 
we decided to start with strongly perturbed cores and examine the 
effect of magnetic fields in such cases. We performed simulations 
of molecular cloud cores with initial m = 2 density perturbations 
with amplitudes of 10, 20 and 50 percent. Initial conditions with 
50 percent density perturbations showed essentially no dependence 
of the results on the magnetic field up to mass-to-flux ratios of 
Af /$ « 2 for a field that was initially aligned with the rotation axis 
- all the calculations formed binaries and their separations were 
essentially indistinguishable. Even with a magnetic field initially 
perpendicular to the rotation axis, although there was a decrease 
of the binary's separation with the initial field strength, binaries 
were formed up to mass-to-flux ratios of A//$ ~ 3. This demon- 
strates that if the initial density perturbations are large enough, 
even strong magnetic fields are not able to suppress or even signifi- 
cantly alter the pattern of fragmentation. We note that this result is 
in qualitative agreement with the single MHD calculation of binary 
formation performed by Ziegler (2005). 

As presented in Section 14.21 for spherical molecular cloud 
cores with weaker 10 percent m = 2 density perturbations we 
find that the fragmentation does depend on the field strength and 
the alignment of the field relative to the rotation axis. A stronger 
field results in a closer binary with only a single protostar be- 
ing produced when the the mass-to-flux ratio is reduced to below 
(M/$) « 4 or (Af/$) ^ 8 for initial fields that are aligned or 
perpendicular to the the rotation axis, respectively. However, the 
question arises as to whether this effect is primarily due to the in- 
creased support provided by the magnetic field (i.e. an increase 
in the effective a of the cloud) or due to the loss of angular mo- 
mentum because of the magnetic braking provided by magnetic 
tension forces. It is well established that clouds w ith higher ini- 
tial values of a are less prone to fragmentation (e.g. lMivama et al] 



ll984lTsuribe & Inutsuki 19991) . As we have demonstrated above, 
re-running the m = 2 calculations in which the field is initially 
aligned with the rotation axis without magnetic tension, or indeed 
without magnetic fields at all but with a thermal a (i.e. increased 
sound speed) that is set equal to the effective a in the MHD calcula- 
tions, gives very similar results to using the full MHD equations. In 
particular, we find that the trend of decreasing separation followed 
by the transition to a single protostar as the magnetic field strength 
is increased is almost entirely due to the extra support provided by 
the magnetic field. When the magnetic field and the rotation axis 
are initially perpendicular, the transition to a single protostar occurs 
at lower field strengths than with an aligned field. However, in these 
cases omitting magnetic tension forces results in even weaker field 
strengths producing single protostars. Unexpectedly, magnetic ten- 
sion aids fragmentation in these cases due to magnetic cushioning 
between the two initial density perturbations. In summary, previous 
studies of the effect of magnetic fields on fragmentation have gen- 
erally attributed the suppression of fragmentation to magnetic brak- 
ing via magnetic tension/torsion forces. However, our results show 
that it is not this straightforward. In the calculations presented here, 
magnetic pressure plays the dominant role in suppressing fragmen- 
tation, while magnetic tension either has little effect or aids frag- 
mentation, depending on the field geometry. 



6 SUMMARY 

We have performed magnetohydrodynamic (MHD) simulations of 
the collapse of molecular cloud cores, some with initial m = 2 
density perturbations, and others that were homogenous. 

In terms of the effect of magnetic fields on the fragmenta- 
tion of perturbed molecular clouds, we have two main conclusions. 
First, we find that wide binaries can be readily obtained from per- 
turbed molecular cloud cores even with mass-to-flux ratios as low 
as A//$ = 2 — 3. Since most molecular cloud cores are observed 
to be strongly aspherical, we conclude that magnetic fields may not 
be as significant a problem to binary formation as many past stud- 
ies have suggested. Second, in agreement with past studies, we find 
that magnetic fields act to suppress fragmentation. However, con- 
trary to past studies that have emphasized the importance of mag- 
netic tension forces and magnetic braking in suppressing fragmen- 
tation, we find that the extra support given by magnetic pressure 
over the thermal support is the dominant reason for the suppres- 
sion of fragmentation and that magnetic tension can actually aid 
fragmentation. 

For all calculations, we find that stronger magnetic fields re- 
sult in longer delays to the collapse due to the increased support 
provided by the field above that provided by thermal pressure alone. 
The delaying of the collapse has a potentially crucial effect on pro- 
tostellar discs that form around protostars. Because the infall rate of 
mass onto the disc from the envelope is reduced, the discs are less 
prone to gravitational instabilities. This may have at least three im- 
portant effects. First, high infall rates on to discs have been shown 
in the past to be responsible for driving gravitational instabilities 
that are strong enough to result in disc fragmentation and the pro- 
duction of binary and multiple systems. Thus, we conclude that 
magnetic fields tend to inhibit disc fragmentation. Second, even for 
initial conditions that in the absence of magnetic fields give infall 
rates insufficient to cause disc fragmentation, the discs frequently 
generate spiral density waves that efficiently transport angular mo- 
mentum. With magnetic fields, the lowering of the infall rate on to 
the disc decreases the strength and importance of these instabili- 
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ties. Third, spiral density waves have been suggested as a possible 
mechanism for producing high enough concentrations of metre- 
sized planetesimals to produce kilometre-sized planetesimals via 
gravitational instability. Since we find that magnetic fields decrease 
the strength of the spiral density waves in protostellar discs, mag- 
netic fields may also decrease the likelihood of forming planetary 
systems by such a mechanism. 
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